Finite driving rates in interface models of Barkhausen noise 
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We consider a single-interface model for the description of Barkhausen noise in soft ferromagnetic 
materials. Previously, the model had been used only in the adiabatic regime of infinitely slow field 
ramping. We introduce finite driving rates and analyze the scaling of event sizes and durations for 
different regimes of the driving rate. Coexistence of intermittency, with non-trivial scaling laws, 
and finite-velocity interface motion is observed for high enough driving rates. Power spectra show 
a decay ~ u) , with t < 2 for finite driving rates, revealing the influence of the internal structure 
of avalanches. 
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I. INTRODUCTION 

The Barkhausen effect p] constitutes a useful, non- 
destructive probe into the domain structure of soft fer- 
romagnetic materials. By ramping an externally applied 
magnetic field within appropriately chosen ranges of in- 
tensity and driving rate, one causes microscopic domain 
wall motions, i.e. magnetization jumps, within a sample. 
The consequent changes in magnetic flux, in turn, in- 
duce a time-dependent electromotive force V(t) on a coil 
wrapped around the sample. Analysis of V(t), assisted by 
suitable theoretical modeling, may provide insight both 
into the domain structure itself and its dynamical behav- 
ior. 

Recent increased interest in the problem stems partly 
from the successful application of methods of non- 
equilibrium statistical mechanics, such as a Langevin de- 
scription via Fokker-Planck equations and avalanche 
models p-^4[. Though there is general agreement on 
the basic underlying mechanisms of Barkhausen noise, 
pinning down specific details to features of theoretical 
models has proved a complex task. For instance, while 
several formulations have focused on the motion of a sin- 
gle interface in a disordered medium |g,||,gj|] , others have 
adopted a picture of nucleation of multiple domains in a 
random- field Ising system Q. 

The rapid and discrete variation of the magnetization 
observed in experiments is a direct manifestation of the 
existence of different time scales. Often the magnetiza- 
tion jumps are regarded as instantaneous events, but this 
is a simplification that allows the description of a limited 
dynamical regime, as finite driving rates are known to 
affect several aspects of the hysteresis cycle. 

Barkhausen noise is observed in the central part of the 
hysteresis cycle, near the coercive field where the magne- 
tization process is mainly due to domain wall motion, so 
interface models are usually successful in describing the 
associated non-trivial scaling laws pMpL[L2J14|. In this 
paper we study the evolution of a single two-dimensional 
interface as it advances driven at a finite rate. The scal- 



ing of event sizes and durations is analyzed for different 
regimes of the driving rate. 



II. MODEL; FINITE DRIVING RATES 

Here we shall use the single-interface model introduced 
in Ref. j| for the description of Barkhausen noise, with 
adaptations for a finite driving rate. The interface at 
time t is described, in space dimensionality d, by its 
height h(pi,t), where pi is the position- vector of site i 
in a (d — l)-dimensional lattice. Here we consider only 
d = 3. At each t, the height function hi = h(p~i,i) 
is assumed to be single- valued, thus the interface ele- 
ment corresponding to the d-dimensional position- vector 
fl = (pi, hi) may be unambiguously labeled by i. Simu- 
lations are performed on an L d ^ 1 x oo geometry, with the 
interface motion set along the infinite direction. There- 
fore finite-size effects are controlled by the length param- 
eter L M . Each element i of the interface experiences a 
force of the form: 



/, = u(fi) + 



3=1 



where 



H e = H -ijM . 



He , (1) 



(2) 



The first term on the RHS of (nf) represents the pinning 
force, u, and brings quenched disorder into the model by 
being chosen randomly, for each lattice site ft, from a 
Gaussian distribution of zero mean and standard devia- 
tion R. Large negative values of u lead to local elements 
where the interface will tend to be pinned, as described 
in the simulation procedure below. The second term cor- 
responds to a cooperative interaction among interface el- 
ements, assumed here to be of elastic (surface tension) 
type. In this term, £ 3 (i) is the position of the j-th near- 
est neighbor of site i and 7 is the coordination number of 



the (d — l)-dimensional lattice over which the interface 
projects. The tendency of this term is to minimize height 
differences among interface sites: higher (lower) interface 
elements experience a negative (positive) force from their 
neighboring elements. The force constant k gives the in- 
tensity of the elastic coupling, and is taken here as the 
unit for /. The last term is the effective driving force, 
resulting from the applied uniform external field H and 
a demagnetizing field which is taken to be proportional 

rd—1 

to M = (1/L ) Yli=i hi, the magnetization (per site) 
of the previously flipped spins for a lattice of width L. 

For actual magnetic samples, the demagnetizing field 
is not necessarily uniform along the sample, as implied 
in the above expression; even when it is (e.g. for a uni- 
formly magnetized ellipsoid) , r\ would depend on the sys- 
tem's aspect ratio. Therefore, our approach amounts to a 
simplification, which is nevertheless expected to capture 
the essential aspects of the problem. See Ref. [O] for a 
detailed discussion. 

Here we use R = 5.0, k = 1, r\ = 0.05, the same values 
as in the d = 3 simulations of Ref. || . 

We start the simulation with a flat wall. All spins 
above it are unflipped. The applied field H is set to the 
saturation value of the effective field H e , in order to min- 
imize transient effects (see e.g. Figure 1 of Ref. ||). The 
saturation H e depends on R, k and r\ {not noticeably on 
L) || , and can be found from small-lattice simulations. 
The force fi is then calculated for each unflipped site 
along the interface, and each spin at a site with fi > 
flips, causing the interface to move up one step. The 
magnetization is updated, and this process continues, 
with as many sweeps of the whole lattice as necessary, 
until fi < for all sites, when the interface comes to a 
halt. The external field is then increased by the mini- 
mum amount needed to bring the most weakly pinned 
element to motion. The avalanche size corresponds to 
the number of spins flipped between two interface stops. 

In line with standard practice [[7|-pd| our basic time 
unit is one lattice sweep, during which the external field 
is kept constant, and all spins on the interface are probed 
sequentially as described above. In the adiabatic regime, 
the external field is kept constant for the whole duration 
of an avalanche, i.e. for as many sweeps as it takes un- 
til no unstable sites are found along the interface. It is 
then increased by the amount needed to flip the weakest 
one. At finite driving rates, the field is increased by a 
fixed amount, henceforth denoted A, at the start of each 
sweep while an avalanche is taking place. Eventually, no 
more unstable sites will be left, and then one proceeds 
as in the adiabatic regime, increasing the field as much 
as necessary to start a new avalanche. Calling SH this 
latter quantity, the time interval between the end of one 
event and the start of the next is then SH/A. 

In what follows, we usually collected samples of 10,000 
avalanches for each simulation, and data analysis has al- 
ways been performed with the whole sets of data; how- 
ever, in figures such as scatter plots of duration versus 



size we display only a representative subset, typically 
500-1,000 events, in order to avoid unnecessary clutter. 



III. SIZE AND DURATION DISTRIBUTIONS 

We consider the single-interface model introduced in 
Ref. p|, initially in the adiabatic limit. In Ref. || we 
showed that, in the context of this model, the upper cut- 
off in the avalanche size distribution is simply a mani- 
festation of finite-size effects, i.e. of the width L of the 
cross-section of our simulated systems perpendicular to 
the direction of interface advance. Fitting our data to the 
customary power-law form multiplied by an exponential 
decay, namely 



P(s) oc s T exp(— s/sq) 



(3) 



we had r in the range 1.2 — 1.4 and Sq oc L u , u = 1.4±0.1. 

The above value for u is consistent with l/<7fc ~ 2/3 
[ defined via so ~ fc~ 1 ' <Jfc , k being a generic demagne- 
tization factor], found both from renormalization-group 
analysis and numerical simulations of the interface model 
of Refs. P,0,0|, provided one makes the identification 
k = 1/L 2 , as pointed out in Ref. Q. 

Wc first investigate whether finite widths affect prop- 
erties related to avalanche duration. 
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FIG. 1. Scatter plots of avalanche duration T (number of 
lattice sweeps) versus avalanche size (number of flipped sites), 
in the adiabatic regime. Bottom to top: L = 40, 80, 160, 320; 
500 events displayed for each lattice L. Plots successively 
shifted upwards by a factor of 10 on vertical scale, to avoid 
superposition. Straight line has slope 0.58 (see text). 



In Figure [I one sees that in the adiabatic regime, apart 
from the L-dependent upper cutoff just mentioned, there 
is no distinguishable influence of finite lattice width on 
the distribution of avalanche durations T against size s. 

For small events with s < 10 flipping sites, dura- 
tion varies approximately linearly with size (collections 
of independent, localized flippings). For larger, collec- 
tive events, non-trivial scaling takes hold; the relation- 
ship is described by the power law T ~ s avz (notation 
borrowed from Refs. ||,[fl]]). Least-squares fits of data 
excluding small avalanches up to s m i n = 10 — 20 give 
avz = 0.58 ± 0.01, very similar to that previously ob- 
tained in simulations of the nucleation model, variously 
quoted as 0.57±0.03 (Ref. §) or 0.58 (Ref. 0). Thus, in 
this aspect at least, there is universality between single- 
interface and nucleation pictures. 

Similar values of avz are obtained also from slightly 
different variants of interface models fl5],|l2]-[l5| . In partic- 
ular, simulations reported in Ref. |T|] give z — 1.56±0.06 
and the roughness exponent ( = 0.75 ± 0.02, which, to- 
gether with the scaling relation a^ 1 = v{2 + £) (see Eqs. 
(34)-(37) of Ref. @), yields avz = 0.57 ± 0.02 . 

Furthermore, the distribution of durations fits reason- 
ably well to a power law with exponential tail, similarly 
to the size distribution Eq. (j|): 



P{T) <x T- a exp (-T/T ) , 



(4) 



where, from standard probability theory, a = 1 + (r — 
l)/avz (— 1.5 ± 0.2 from the values of r and avz quoted 
above). Indeed, analyzing the data shown in Fig. Hone 
gets a in the range 1.3 — 1.5. The cutoff T is expected 
to scale with the size cutoff sq a s Tq ~ s 



q vz , therefore 
(with s ~ L 1A±{)1 |6|), one must have T ~ £0.8i±0.06 
Direct analysis of data gives the latter exponent varying 
in the range 0.65 — 0.9, broadly compatible with this. 

Next, we gauge the effects of varying driving rates on 
the size-duration relationship. For sufficiently fast driv- 
ing rates, one expects coalescence of avalanches which 
would be separate events in the adiabatic regime. 

We start by fixing L — 80. A, the external field 
increase at the start of each lattice sweep during an 
avalanche, is the driving rate . In Figure @ one has 
A = 10~ 5 , 5 x 10~ 5 , 7.5 x 10~ 5 , and 10" 4 . 

The plot for A = 10~ D is identical, apart from small 
fluctuations, to the corresponding one for L — 80 in the 
adiabatic regime, shown in Fig. |]. We have checked that 
the same happens for the in-between driving rates 10 _m , 
to = 8, 7, 6 (as expected). Again, least-squares fits give 
avz = 0.58 ±0.01 not only for the whole set of A = 10~ D 
data, but also for the initial parts (with s < 10 5 ) of 
those with faster driving rates. These latter plots will be 
discussed further ahead. 

We conclude that, for L = 80 and within the range of 
A = — 10~ 5 , there is no influence of the driving rate on 
the T — s relationship, including the maximum avalanche 
size, which remains at s max ~ 10°. So far, the results: 
(i) t = 1.3 ±0.1; (ii) a = 1.5 ±0.2 and (iii) independence 



of behavior on driving rate (at least, within fairly well 
defined windows of A, for given L) show that the model 
under consideration shares the same universality class of 
the interface model discussed in Refs. P,|l2|-[l4| , when 
dipolar interactions are neglected. 
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FIG. 2. Scatter plots of avalanche duration T versus 
avalanche size, for L = 80 and different driving rates A. Bot- 
tom to top: A = 10" 5 , 5 x 10" 5 , 7.5 x 1(T 5 and 10" 4 ; 500 
events displayed for first three, 1, 000 for A = 10 -4 . Plots suc- 
cessively shifted upwards by a factor of 10 on vertical scale, to 
avoid superposition. Full straight line has slope 0.58; dashed 
line has unitary slope (see text). 



IV. RATE-DEPENDENT BEHAVIOR 

We now investigate the different picture appearing for 
higher driving rates, already shown in Fig. @. The scatter 
plot of duration against size for A = 5 x 10~ 5 shows 
traces of a different behavior for s > s max as a higher 
slope develops at that point. As A is further increased, 
it can be seen that there is coexistence between non- 
trivial scaling for sizes < the adiabatic cutoff, and T ~ s 
behavior (i.e. interface motion with a finite velocity) for 
larger avalanche sizes. 

A qualitative explanation of the above results goes as 
follows. With Hf(i) and Hi(i + 1) being, respectively, 
external field at the end of the i-th avalanche and at the 
start of the i + 1-th, typically the corresponding gaps 
5Hf d = Hi(i + 1) — Hf(i) in the adiabatic regime are 
larger than, say, several times A for A < 10 -5 ; therefore 
few avalanches merge for such driving rates. For larger 
A more gaps are closed, and the distribution changes 
significantly. 

In order to gain a quantitative understanding of this, 
at the same time checking for a possible L-dependence, 



we study the probability distribution of SH ad for different 
lattice widths. 

In Fig. we show the cumulative probability 
P(SH ad < 6H Q ) of 5H ad being smaller than 8H . Before 
analyzing the L— dependence of the curves, we focus on 
L = 80. One sees that P(SH ad < SH ) - 3 x 1CT 3 ,7 x 
1(T 2 , respectively for SH = 1CT 5 , 1(T 4 . From Fig. §, the 
maximum avalanche duration is T max ~ 5 x 10 3 , suggest- 
ing that the finite driving rate is irrelevant as long as the 
quantity V = T(L,A) = T max P{6H ad < A) < O(10). 
On the contrary, for A = 10~ 4 Figs, p] and H show that 
V ~ 0(100). 

The inset of Fig. | shows that P(SH ad < 5H ) is, 
to a good approximation, a function of L x SHq, where 
the best collapse plots are obtained in the range x = 
1.9 ±0.1. Though at this point we are not able to advance 
an argument, it may be that x — 2 exactly. 
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FIG. 3. Adiabatic regime: probability P(SH ad < SH ) of 
the gap between consecutive avalanches being smaller than 
SHo, against SH . Left to right: L = 320, 160, 80 and 40. 
40,000 events for L = 40, 10,000 for larger widths. Inset: 
P(5H ad < SH ) against (L/40)" 8H , x = 1.9. Same ranges 
as main plot, on both horizontal and vertical axes (see text). 



concern in the study of equilibrium second-order phase 
transitions, of trying to extract the thermodynamic limit 
from finite-size scaling: finite- L results must be analyzed 
in their own right. Nevertheless, it is crucial to investi- 
gate the size dependence of relevant quantities, as done 
here, in order to separate truly L-independent features 
from crossover behavior. 

The departures from rate- independent behavior, as de- 
picted in Fig. ||, require further analysis. A least-squares 
fit of data in that Figure for A = 10~ 4 , 10 < s < 5 x 10 4 
(the horizontal extent of the full line shown there) gives 
avz = 0.58 ± 0.01. The events in this range of s are 
~ 20% of the total number of avalanches; another 40% 
are small ones with s < 10, with the remaining 40% 
larger than 5 x 10 4 sites. The overall probability density 
P(s) for A = 10 -4 , with the customary logarithmic bin- 
ning, is displayed in Fig. |4J, for the four driving rates of 
Fig. H . The full straight line suggests that the A = 10~ 4 
data with 1 < s < 5 x 10 can be fitted by a power law 
P(s) oc s~ 16 . This larger effective value of r can be un- 
derstood by observing the depletion of events with sizes 
10 4 — 10 5 in the cross-over region of Fig. 0. Curves for 
A = 5 x 10 -5 and 7.5 x 10 -5 display intermediate behav- 
ior, with only an incipient shoulder at 10 < s < 10 5 for 
the former. 
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It is thus clear that as L increases the range of A for 
which adiabatic behavior dominates is shrunk. Recall 
that, on top of the scaling of probabilities with L x SHq, 
T m ax (which scales with the cutoff To) also depends on 
L v , v ~ 0.81, see paragraph below Eq. ([|). Thus, in the 
limit L — > oo one should have adiabatic-like properties 
only strictly at A = 0. Before attaching much signif- 
icance to this, one must bear in mind the main result 
of Ref. § : in the present model, L is closely correlated 
with the (average, or typical) domain size in experimen- 
tal situations. Therefore, one does not have the usual 



FIG. 4. Probability distributions of avalanche sizes for 
L = 80. Bottom to top: A = 10~ 5 , 5 x 10~ 5 , 7.5 x 10~ 5 and 
10 -4 . Curves successively shifted upwards by a factor of 10 
on vertical scale, to avoid superposition (baseline is A = 10~° 
plot). Straight lines have slopes as indicated, and are intended 
as guides to the eye. 

We note that in earlier studies, both experimental |lq ] 
and numerical p7[ , it was found that the exponents r of 
Eq. (Q) and a of Eq. ([|) decrease as driving rate increases, 
apparently in contradiction with the above. Also, a small 



shoulder can be seen in the data for alloys under stress 
in Figure 1 (a) of Ref. |14|] . We comment on each of these 
in turn. 

Firstly, recent work |L4| shows that the FeSi alloys 
of Ref. p6[ and, e.g., amorphous alloys under stress 
belong to different universality classes; while the for- 
mer follow a mean-field description (giving rise to ac- 
tual rate-dependent exponents), the latter display rate- 
independent exponents with values close to those found 
both for the present model, and for the model of 
Refs. |||l2] H when dipolar interactions are neglected 
(see remarks at the end of Sec. Ill) . 



As regards the simulations of Ref. |17|, it may be that 
the model used here does not belong to the same dynam- 
ical universality class as the (two-dimensional) random- 
field model with vacancies used there. However, exami- 
nation of the data in Fig. |j suggests an alternative expla- 
nation, as follows. Though fits of the straight-line parts 
of the distributions give r increasing with A, fits of the 
whole sets of data to a form such as Eq. p) , with two 
free parameters (r and s ) i n f ac t gi ye decreasing values 
of r for larger A, on account of the large- s shoulders . 
Of course, this happens at the expense of the quality of 
the fit; however, for less dramatic departures than those 
shown in the Figure, the loss of quality might not be 
obvious. 

Finally, the incipient shoulders shown in Ref. [M ap- 
pear to be bona fide candidates for the above descrip- 
tion, as they correspond to materials to which, so far, 
the present model seems to fit well. We believe that a re- 
examination of experimental and simulational raw data, 
in search of a coexistence of regimes, would be worth- 
while. As pointed out in Ref. |L7|, "[a faster driving rate] 
overdrives weaker pinning centers thus rendering the oc- 
currence of larger avalanches more probable". For the 
present model this indeed happens, only it does so at 
the expense of depleting the histograms of occurrence 
of events < the respective maximum size for the adia- 
batic regime. We expect the present study to motivate 
further experimental and numerical investigation along 
these lines. 



tion will lead to y — 1 . Figure || shows the plot of 1 — y 
versus A for L — 80. The best non- linear fit of the whole 
set of points is given by 1 — y = cxp[— (A/Ao) a ] with 
the optimum values of the free parameters: a — 1.2 and 
Aq = 1.57 x 10~ 5 . The inset shows the rate-independent 
regime corresponding to y <C 1, where one clearly has 
y ~ A. Such linearity is to be expected, as essentially the 
same events occur for any A in this interval: for a given 
avalanche, Hf(i) — Hi(i) is proportional to A, while the 
denominator Hf(N) — Hi(l) is unchanged. The fact that 
the best- fitting value of a is 1.2 indicates that the decay of 
1— y for large A is in fact faster than exp[— (A/Ao) 1 ' 2 ], so 
one is having a compensation among the different regions 
of the plot, in order to minimize the overall deviation. A 
fit of the subset of data for A > 10 -5 indeed gives a ~ 2. 




FIG. 5. Plot of 1 — y (y is the fraction of the total du- 
ration spent during avalanches) against driving rate A, for 
L = 80. Solid curve is fit to 1 - y = exp[-(A/A ) a ], a = 1.2, 
Ao = 1.57 x 10 . Inset: y against A for slow driving; dashed 
line has unitary slope. 



V. INTERMITTENCY 

It is also important to analyze how the intermittency of 
events is gradually lost as more and more avalanches co- 
alesce with increasing driving rate. This can be done by 
defining y as the fraction of time spent during avalanches. 
The duration of an avalanche being given by T(i) — 
(Hf(i) — Hi(i))/A, the overall duration of a simulation 
with N events is T N = (H f (N) - H t (l))/A, therefore 

y = (i/T N )EliT(i) = Eti(tf/(*)-tfi(*)]/[#/W- 

Hi(l)]. In the adiabatic regime A = 0, avalanches are in- 
stantaneous, and y — 0. As A is increased, avalanches 
will be observed part of the time, so < y < 1, and in 
the limit of large A, one expects that a depinning transi- 



VI. POWER SPECTRA 

Power spectra and their analysis are an important tool 
for the understanding of noise in disordered systems; 
specifically for the Barkhausen effect, see e.g. Ref. |ll| for 
a recent survey of results and references; see also Ref. [H . 

We shall always make L = 80 in this section. With the 
unit of time given by a lattice sweep, the maximum fre- 
quency to be analyzed is u>m — 1/2, as standard Fourier 
theory prescribes. 

We attempt to concentrate on the non-trivial scaling 
regime. To this end, we recall from Fig. |]that, whatever 
the driving rate, events of size s < 10 sites (of which there 



is always a large fraction of the total) are collections of 
independent, localized dippings. Thus their duration is 
proportional to size. Therefore, when looking at power- 
spectrum data with frequencies u> ^ 0.1, one will have a 
strong input from such non-critical flippings. Though one 
might think of ways to expurgate the respective contribu- 
tions, we shall simply restrict ourselves to the frequency 
region w < 0.1. As shown in the following, this still 
leaves a suitably wide window of observation, in what we 
call intermediate-frequency range (for u) — > the power 
spectrum goes flat, as details of the temporal series are 
washed out on long time scales). 

The discussion of power spectra revolves mainly 
around the power t with which decay sets in, in the 
intermediate-frequency range: P{oj) ~ w . The follow- 
ing points were made in Ref. pT| : (1) While one has £ = 2 
for white noise and mean-field descriptions, and early 
studies of the Barkhausen effect predicted t = (3—r)/avz 
(= 2.9 ± 0.2 for the present model, with r = 1.3 ± 0.1, 
avz = 0.58 ± 0.01), it was found that for r < 2 one 
should have t — 1/avz {— 1.72 ± 0.03 here), instead 
of (3 — t)I avz. Several experimental and simulational 
results were shown to be compatible with the latter find- 
ing. (2) It was remarked that the result of earlier cal- 
culations of the power spectrum |ll| for sandpile mod- 
els, which gave P(ui) ~ u>~' 2 and were cited as an ex- 
planation for such behavior in (among other systems) 
Barkhausen noise, depends crucially on the assumption 
that the avalanche shape can be approximated by a box 
function: V(t) = S/T (0 < t < T) for an avalanche of 
size S and duration T. 

In the present model, one can tune the degree to which 
the internal structure of an avalanche is taken into ac- 
count, by varying A. While for A = all events are seen 
as spikes of zero duration, the internal fluctuations within 
avalanches become much more noticeable as A increases, 
even still within the adiabatic regime. Recall that, for 
A = 10~ 5 , L = 80, avalanches take up ~ 40% of the 
total duration of a simulation, see inset of Fig. g. 

Below, we set out to probe points (1) and (2). Accord- 
ingly, in Fig. |6| we plot the power spectra for A = 10~ 7 
(deep within the adiabatic regime), 10~ 6 , 10~ 5 and 10~ 4 . 

Starting from the slowest driving rate considered, one 
sees that disregarding the internal structure of avalanches 
indeed yields a dependence of the power spectrum on 
ui~ 2 . This is entirely consistent with point (2) above. The 
next graph, A = 10 -6 , is somewhat difficult to interpret 
on its own. However, the trend becomes clearer when the 
A = 10~ 5 data are taken into account: as more details of 
the intra-avalanche fluctuations enter into the spectrum, 
its decay becomes slower, ~ ll>~ 1,5 . Though the numeri- 
cal values are somewhat off the mark, one definitely sees 
that the trend is towards t = I /avz when intra-avalanche 
correlations are considered, as opposed to the alternative 
t = (3 — t)/ avz. This is in support of point (1) above, 
showing that very likely the present model behaves simi- 
larly to the random- field one of Ref. flit], in this respect. 

Finally, for A = 10 there is an apparent trend re- 



versal; from a least-squares fit of data in the range of u 
shown, one has the exponent 1.83 ± 0.01. It must be re- 
called that one is then clearly in the coexistence regime 
explained above (that is, away from a purely intermit- 
tent, adiabatic framework); therefore there is strong in- 
fluence of inter-avalanche correlations. It is then not sur- 
prising that the picture starts to differ e.g. from that 
of Ref. [[Ll| , where only inter-avalanche correlations were 
taken into account. 

Clearly, more work is needed to sort out this latter 
point. As one goes deeper into the depinned regime, it 
may well be that the power spectrum decay returns to the 
lo~ 2 form characteristic of uncorrelated events. In such a 
scenario, the above value of t = 1.83 would in fact be an 
effective exponent, marking a crossover towards t = 2. 
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FIG. 6. Continuous lines: power spectra for L = 80 and 
various driving rates. Bottom to top, A = 1CP 7 10~ 6 , 1CP 5 
and 10~ 4 . The vertical coordinates have been shifted by dif- 
fering amounts for each plot, so they could fit in a single 
graph. Dashed lines have slopes as indicated, and are in- 
tended as guides to the eye. 



VII. CONCLUSIONS 

We have used a single-interface model with an ad- 
justable finite driving rate to simulate the time sequence 
of Barkhausen jumps of a three-dimensional system. Di- 
rect observation of duration as a function of size deter- 
mines the existence of two dynamical regimes regard- 
ing the avz exponent, that is, for low driving rates 
(A < A C (L)), avz — 0.58 up to a limiting event size, 
and the avalanche dynamics is basically rate indepen- 
dent. For higher driving rates (A > A C (L)) the previous 
regime coexists with one for which avz = 1. The rate 



dependency of the second regime is evident as we ana- 
lyze the probability distribution of avalanches sizes: for 
high driving rates it deviates considerably from the usual 
form, P(s) oc s~ T exp(— s/sq). The passage from one 
regime to the other is rather sharp, and the correspond- 
ing value of driving rate, A c , depends on the system size: 
an infinite system will have a rate-independent dynamics 
for A = only, that is, A c (oo) = 0. 

Considering only the power-law portion of the P(s) 
graphs for finite driving rate, we find exponents that in- 
crease with the driving rate, in apparent contradiction 
with previous theoretical and experimental results. On 
the other hand, when the fitting assumes a power law 
with an exponential cut-off, and the whole set of data 
is taken into account, the effective value of r decreases 
as the driving rate is increased. It is clear though that 
this form does not provide an adequate description of the 
simulated data for A > A c , so we believe that this may 
be at least part of the explanation for inconsistency in 
previously reported values of r. 

The power spectra for various driving rates clearly 
show that with increasing driving rates, intra-avalanche 
correlations become more relevant, as the time scale in- 
volved reveals details of events occuring with a finite du- 
ration. A direct consequence is the relation P(lo) ~ w~* 
observed: as the internal structure of the avalanches is 
probed, t decreases from the value t — 2, characteristic 
of adiabatic time series. 

In summary, we have studied the effect of a finite driv- 
ing rate in the scaling properties of the Barkhausen noise. 
As our ultimate goal is the description of experimental 
results, it is important to understand the limitations in- 
volved in real experiments as compared to computational 
ones. In principle simulations may use any value for the 
driving rate, as well as any system size, or at least we 
may say that our choice of values is broad as compared 
to real experiments. A typical experiment usually has 
driving rates spanning only one order of magnitude, and 
values of domain sizes predetermined by the fabrication 
conditions of the sample and applied stress. So, as the 
experiment is designed with these parameters, the scaling 
regime is basically already chosen. Any further compari- 
son with theoretical or simulation results must be careful 
in the sense that the same regime has to be studied. 
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